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Knowledge of atmospheric muon flux intensity at large depths is extremely important for neutrino 
telescopes located deep under ground, water or ice. One of the methods to transform muon sea- 
level spectrum into depth one is to apply Monte Carlo technique which directly takes into account 
stochastical nature of energy loss. In order to decrease computation time down to acceptable level 
one has to use simplifications resulting in systematic errors which in some cases may distort result 
essentially. Here in this paper we present our analysis for dependence of computed depth muon flux 
upon the most important parameters of muon transport Monte Carlo algorithm which was done with 
the MUM (MUons+Medium) code. Contribution of different simplifications to the resulting error 
is considered, ranked and compared with uncertainties which come from parametrization accuracy 
both for sea-level muon spectrum and for muon cross sections. 

PACS number(s): 13.85Tp, 96.40.Tv, 02.70Lq 



I. INTRODUCTION 



In order to obtain muon spectrum and value for integral muon flux at large depths where existing and planned 
neutrino telescopes are located (one to several kilometers of water equivalent) one has to perform several steps starting 
with interactions of primary cosmic rays in atmosphere. The last step of this procedure is transportation of muons from 
sea level down to detector location. It can be done both with analytical or semianalytical technique and Monte Carlo 
(MC) simulation which directly takes into account that muon energy loss represents an essentially stochastic process. 
MC codes for muon transportation have to use some simplifications to optimize equilibrium between reasonable 
computation time, desirable accuracy and neccessary statistics. 

The usual approximation when calculating the muon propagation through a medium is the following: muon inter- 
actions with comparatively large energy transfers, i.e., when fraction of energy lost v — AE/E exceeds some value 
v cut ~ 10~ 3 -10~ 2 , are taken into account by direct simulation of v > v cu t for each single interaction according to 
shape of differential cross sections da(E,v)/dv (these interactions lead to "stochastic" energy loss or SEL) while the 
part of interaction with relatively small v is treated by the approximate concept of "continuum" energy loss (CEL), 
i.e., using the function [dE(E)/dx]cEL which is obtained from formula 
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Here Na is the Avogadro number; A is the atomic weight; indexes 7, e + e~, pn and S correspond to bremsstrahlung, 
direct e + e~-pair production, photonuclear interaction and knock-on electron production, respectively; v m i n and v max 
are the minimum and the maximum kinematically allowed fraction of energy lost for corresponding process; term 
[dE{E)/dx\B-B represents Bethe-Bloch formula with correction for density effect. Notice that actually one has to 
compute CEL, as well as all total cross sections separately for each kind of atoms given material consists of and then 
add them to each other with different weights but here and below in the text we omit this detail and give only general 
expressions. One is forced to decompose energy loss because simulation of all interactions with v > v min would result 
in infinite computation time due to steep dependence of muon cross sections on v (they decrease with v at least as 
do~(E,v)/dv cx u _1 and for some processes are not finite at v — > 0). Number of interactions to be simulated per unit 
of muon path grows, roughly, as iV, n t cx v~ ut along with computation time. On the other hand, setting v c ut to large 
value may affect simulation accuracy. Thus, the question is how large value ofv cu t may be chosen to keep result within 
desirable accuracy? 
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The second problem can be formulated as follows: if it is necessary to include ionization in SEL, at all? Small 
energy transfers strongly dominate at knock-on electron production (da s (E , v) / dv oc v~ 2 ), so this pro cess is almost 
non-stochastic and it seems to be reasonable to omit the last term in expression for CEL (Eq. (1.1)) and exclude 
knock-on electrons from simulation procedure when simulating SEL, i.e., to treat ionization completely as "continuous" 
process which saves computation time noticeable. How much does it affect the result of simulation? 

Influence of these factors on simulated result was discussed in literature (see, e.g., Refs. fl|-|]]) but, in our opinion, 
more detailed analysis on which further improvement of muon propagation MC algorithms could be based is still 
lacking. In presented work we undertook an attempt to investigate in details with the MUM (MUons+Medium) muon 
propagation code (Ref. ||) the influence of v cut and model of ionization energy loss upon result of MC simulation. 
Resulting errors are compared with uncertainties which come from parametrization accuracy both for sea-level muon 
spectrum and for muon cross sections. Our analysis allows to rank different uncertainties according to their importance, 
expose ones to which special attention should be put and choose the most adequate setting of MC parameters at 
simulation for this or that purpose. 

We briefly describe the main features of the MUM code which was a basic instrument for reported investigation 
and give a short review of principal variables have been studied in Sec. [o|. Sec. ^ presents results of simulation for 
survival probabilities and muon flux intensities with different settings of simulation paremeters. In Sec. ^ we analyse 
results of simulations presented in Sec. Ill, compare them with other published results and give some basic inferences. 
Sec. |v| gives general conclusions. 



II. METHOD 

Description of the MUM code which was used to obtain presented results has been published in Ref. [^) . Here we 
only dwell briefly on those features of the code which seem to be important in the frame of discussed problem: 

- the most recent results on parametrizations for muon cross sections are used in the code; 

- we tried to decrease the "methodical" part of systematic error due to interpolation, numerical integration, etc., 
down to as low level as possible, especial attention was put on simulation procedures for free path and fraction 
of energy lost; 

- the most important parameters, like value of v cu t, model for ionization loss, kind of medium, parametrizations 
for muon cross sections are changeable and represent input parameters for initiation routine; 

The accuracy of simulation algorithm used in MUM was shown in Ref. H to be high enough to perform systemat- 
ically significant analysis which is presented below. 

In order to study how different factors influence upon result of simulation we performed several sets of simulations 
both for propagation of monoenergetic muon beam and muons sampled by real sea-level spectrum (in the later case 
we limited ourselves by simulation only vertical muons) through pure water down to depths from D = 1 km to D = 
40 km. Of course, depths of more than several kilometers of water for vertical muons do not concern any real detector 
but simulations for larger depths allow us to study general appropriatenesses for large muon ranges which correspond, 
for instance, to nearly horizontal muons. Several runs were done for standard rock (A = 22, Z = 11, p = 2.65 g 
cm' 3 ), as well. Muons whose energy decreased down to E = 0.16 GeV (the Cherenkov threshold for muon in water) 
were considered as stopped ones. We tested different settings of parameters which were as follows. 

(a) v cu t, which changed within a range of v cut = 10~ 4 to v cut — 0.2. Actually, "inner" accuracy of the MUM 
code becomes somewhat worse at v cut > 5xl0~ 2 , especially if fluctuations in ionization are not simulated 
(Ref. §). So, results for Vcut = 0.1 for v cut = 0.2 are presented here only to illustrate some general qualitative 
appropriatenesses. 

(b) Model for ionization loss which was treated both as completely "continuous" and "stochastic" with simulation 
of energy lost if knock-on electron energy AE > v cut E. 

(c) Parametrization for vertical sea-level conventional atmospheric muon spectrum. Two spectra were tested, namely 
one proposed in Ref. H (basic): 
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and widely used Gaisser spectrum (Ref. Q): 
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(d) Parametrization for total cross section for absorbtion of a real photon of energy v = s/2tojv = vE by a nucleon 
at photonuclcar interaction which was treated both according to Bezrukov-Bugaev parametrization proposed in 
Ref g (basic): 

a lN = [114.3 + 1.647 ln 2 (0.0213 v)] fib (2.3) 

and ZEUS parametrization (Ref. ||): 

a lN = ( 63.5 s 097 + 145 s~ a5 ) fib (2.4) 

(y and s in GeV and GeV 2 , correspondingly). 

(e) A factor k a which all muon cross sections were multiplied by to test influence of uncertainties in cross sections 
parametrization (and, concequently, in energy loss) upon result. We used k a — 1.0 as a basic value but in some 
cases set also k a = 0.99 and k a — 1.01, which corresponds to decrease and increase of total energy loss by 1%, 
respectively. Note that it is an "optimistic" evaluation, the real accuracy of existing parametrization for muon 
cross sections is worse (see Refs. [ |10|]ri| ). 

For each run we fixed the muon spectra at final and several interim depths. The differences between obtained 
spectra were a point of investigation. 

III. RESULTS 

A. Propagation of monoenergetic muon beams 

At the first set of simulations we propagated monoenergetic muon beams of 4 fixed initial energies E s = 1 TeV, 10 
TeV, 100 TeV and 10 PeV down to depths D — 3.2 km, 12 km, 23 km and 40 km, respectively, through pure water. 
In each case propagation of 10 6 muons was simulated. Fig. |Tj shows resulting survival probabilities p = No /N s (where 
N s = 10 6 is initial number of muons and Nu is number of muons which have survived after propagation down to 
depth D) vs. v cu t for final and five interim depths. Two cu rves are given on each plot for two models of ionization. 



Also results for k a — 1.00 ± 0.01 and for parametrization (2.4) are presented as simulated with the most accurate 
value v cu t = 10 -4 . 

The following appropriatenesses are visible. 

(a) In most cases with the exception of some plots of the lower row and the left column in Fig. |l| (which corresponds 
to low survival probabilities and low muon initial energies, respectively) uncertainty in our knowlege of muon 
cross sections gives the principal effect which essentially exceeds ones from other tested parameters. 

(b) The difference between survival probabilities for two models of ionization is the less appreciable the larger muon 
energy is. It is quite understandable because at muon energies E < 1 TeV ionization represents the great bulk 
of total energy loss, and vice versa, it becomes minor at E > 1 TeV. Thus, contribution which is given by 
ionization at higher energies is small and, the more, its fluctuations do not play an important role. For muons 
with initial energies B»l TeV fluctuations in ionization become important only at very last part of muon 
path and "are not in time" to produce some noticeable effect. 



(c) Generally, parametrizations (2.3) and (2.4) do not show a noticeable difference in survival probabilities, in most 



cases it is within statistical error or exceeds it only slightly. 

(d) Increase of v cut gives effect of both signs in survival probabilities: function p(v cut ) grows at the beginning of 
muon path and falls at the last part. The same "both-sign" dependencies are observed for ionization model. 
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(e) For v cut < 0.02 - 0.05 there is almost no dependence of survival probability on v cu t with the exception of very 
last part of muon path where survival probability becomes small. Generally, dependence p(v cut ) is the less 
strong the larger initial muon energy is. 

The last item is illustrated complementary by Fig. ^ and Fig. || which show that for all initial energies E s simulated 
survival probability does not depend, in fact, on v cu t until 90% (for E s = 1 TeV) to 99.5% (for E s = 10 PeV) muons 
have been stopped. 

It was shown above what is result of simulations with different models of ionization and values of v cut . It was a 
special point of interest for us to track how does it influence upon behaviour of survival probability. 



E s =1TeV: 



E s =10TeV: 



E s =100TeV: 



E s =10PeV: 
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FIG. 1. Survival probabilities p — Nn/N s (where N B = 10 6 is initial number of muons and No is number of muons which 
have survived after propagation down to depth D in pure water) vs. v cut . Values of p were obtained as a result of MC simulation 
for monoenergetic muon beams with initial energies E s = 1 Tev (1st column of plots), 10 TeV (2nd column), 100 TeV (3rd 
column) and 10 PeV (4th column). Each column contains six plots which correspond to six depths D (which differs for different 
E s ). Closed circles represent survival probabilities which were simulated with ionization energy loss included in SEL along 
with other types of muon interactions. Open circles correspond to computation with completely "continuous" ionization. Two 
horizontal solid lines on each plot show the value for survival probability computed with all muon cross sections multiplied by 
a factor k a = 1.01 (lower line) and k a = 0.99 (upper line) for v cu t = 10~ 4 . Horizontal dotted lines correspond to v cut = 10~ 4 
and cross section for absorbtion of a real photon at photonuclear interaction parametrized according to r ecen t ZEUS data (Ref. 
M, Eq. (2^)) instead of parametrization proposed by L.B. Bezrukov and E.V. Bugaev (Ref. [B], Eq. (^3)) which is basic in 
the MUM code. Note different scales at Y-axis. 
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FIG. 2. Survival probability p vs. depth D down to which muon beam of initial energy E s = 1 TeV (a), 10 TeV (b), 100 TeV 
(c) and 10 PeV (d) propagates through pure water. On each plot 10 lettered curves which correspond to different values of v cut 
are shown. Meaning of letters is as follows: A - v cu t = 10~ 4 , B - v cut = 2xl0~ 4 , C - v cu t = 5xl0~ 4 , D - v cut = 10~ 3 , E - v cu t 
= 2xl0~ : \ F - v cut = 5xl0~ 3 , G - v cut = 10~ 2 , H - v cu t = 2xl0~ 2 , I - v cu t = 5xl0~ 2 , J - v cu t = 10 _1 . This figure displays 
results were obtained by simulation with ionization loss included in SEL. Statistical errors which cause some unsmoothness of 
curves at small p are not shown. Dependence p upon v cu t becomes noticeable only at the last « 1/8 of muon beam path where 
the majority of muons (from 90% to 99.5%, depending upon E s ) has been stopped. 




FIG. 3. Relation p/po vs. p. p is survival probability for muon beam of initial energy E s = 1 TeV (a), 10 TeV (b), 100 TeV 
(c) and 10 PeV (d) at propagation through pure water with ionization included in SEL as simulated for different values of v cu t ; 
po is survival probability simulated under the same conditions for v cu t = 10~ 4 . Difference in p/po becomes noticeable only at 
small values of p, i.e., at the last part of muon beam path. 
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FIG. 4. Muon spectra at different depths D in pure water resulting from 10 6 muons with initial energy E s = 1 TeV obtained 
by simulation with four models. The first three columns represent spectra obtained with ionization included in SEL for v cu t = 
1CF 4 (1st column), 10 -2 (2nd column), and 0.2 (3rd column). 4th column contains spectra obtained for entirely "continuous" 
ionization and v cu t = 10 -4 . On each plot value of survival probability p is indicated (without statistical error which is negligible 
for the most plots and, on desire, can be easily calculated taking into account that p — Nn/W 6 , where Nd is number of muons 
which have survived after propagation down to depth D). 

Fig. U shows how muon spectrum resulting from monoenergetic muon beam with initial energy E s = 1 TeV transforms 
when propagating through pure water down to the depth of 3.2 km. Simulated results for four settings of parameters 
are presented by four columns of plots. The first three columns represent spectra obtained with ionization included in 
SEL for Vcut — 10~ 4 , 10~ 2 and 0.2. 4th column contains spectra simulated with entirely "continuous" ionization and 
v C ut — 10 -4 . The spectra grouped into the first column represent some "ethalon" because they were simulated with 
the most accurate tuning both for v cu t and ionization model. The first three columns demonstrate that "compactness" 
of spectra at the same depth is the higher the more value of v cut is. Put your attention to the right edge of spectra 
which shifts toward low energies when v cu t increases (it is the most noticeably for v cu t — 0.2). The reason is that 
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at any depth energy of the most energetic muons in simulated beam is determined by CEL. These muons due to 
statistical fluctuations did not undergo interactions with v > v cu t and, consequently, lost energy only by CEL which 
increases when v cut increases. That is why the maximum energy in simulated muon beam is lower for large values 
of v cu t- Fraction of muons which did not undergo an "catastrophic" act with v > v cu t till given depth grows with 
increase of v cu t because free path between two sequential interactions with v > v cu t grows approximately as L oc v cut . 
It leads, in particular, to distinctly visible picks in spectra for v cu t — 0.2 consisted just of muons which lost energy 
only by CEL. Also, some deficit of low energy muons appears if one sets v cu t to a large value. In this case left edge 
of spectrum is provided only with muons which interacted with large fraction of energy lost while for smaller v cu t 
an additional fraction of muons comes here. As a result simulated spectrum of initial monoenergetic muons at given 
depth is more narrow if v cu t is large and, on the contrary, more wide if v cu t is small. 

Now it is easy to understand how value of v cut influences on simulated survival probabilities. When simulated muon 
beam goes through medium loosing energy both in CEL and SEL processes, its spectrum is constantly shifting to the 
left (energy decreases). For v cut = 1CP 4 the left part of spectrum reaches E — at a smaller depth comparing with 
larger v cut (because in this case spectrum is wider) and survival probability starts to decrease. At the same depth 
survival probability for v cu t — 10 -2 and v cu t — 0.2 is still equal to 1. Thus, for the first part of path the survival 
probability is always larger for large v cu t. At some depth (which is equal to approximately 2.8 km for considered case) 
compactness of spectra simulated with large v cut starts to play an opposite role. Due to this compactness and higher 
CEL muons stop faster comparing with accurate simulation. So, at the final part of the beam path simulated survival 
probability for large v cut decreases faster comparing with accurate simulation and, for instance, for v cut = 0.02 the 
rest of muon beam which reaches the depth of D = 2.72 km (37% of initial number of muons) completely vanishes 
within the next 30 m of path, while some fraction of muons simulated with v cut — 10 -4 (0.07%) escapes down to the 
depth of D = 3.2 km. 

Qualitatively the same effect leads to the same results if one treats ionization as completely "continuous" energy 
loss. Again, spectra becomes more narrow since fluctuations in ionization do not work and, as a consequence, survival 
probability becomes significantly higher comparing with simulation with accurate treatment of ionization at the 
beginning of muon beam p ath an d falls down essentially faster at the final part of path. 

Results presented in Sec. Ill A show the significant influence which both model of ionization and value of v cut have 
over survival probability for monoenergetic muon beam. But for practical purposes the more important is how this 
factors do work for real atmospheric muons with a power spectrum? 



B. Propagation of muons sampled according to a power sea-level spectrum 

In Fig. |^ we present intensity of vertical atmospheric muon flux / at different depths of pure water D from 1 km to 
20 km vs. v cut as simulated with muons sampled according to sea-level spectrum (2.1). Simulation continued until 
10 4 muons reached given depth. Curves for two models of ionization are shown for each depth along with results for 



k a = 1.00 ± 0.01 at v cu t — 10~ 4 , parametrization (2.4) at v cu t = 10 -4 , sea-level muon spectrum (2.2) at v cu t = 10~ 4 
and all energy loss treated entirely as CEL (for depths D < 5 km only). 

General appropriatenesses for real muon spectrum are qualitatively the same as observed for monoenergetic muon 
beams. 

(a) For all depths at which neutrino telescopes are located it was found to be better to take into account fluctuations 
of energy loss simulated by any model than to treat energy loss as completely "continuous" : muon flux intensity 
computed with non-stochastical model of energy loss is always less comparing with stochastical model, the 
difference reaches 10 % at 3 km w.e. and 40% at 5 km w.e. (at the depth of 20 km of pure water vertical muon 
flux computed with ignorance of fluctuations is only 10 % of simulated flux). 

(b) Like in a case for monoenergetic beams 1%-uncertainty in muon cross sections plays the principal role for 
resulting error in simulated muon depth intensity. This error has a tendency to grow with depth from ±2.5% at 
depth of 1 km w.e. to ±15% at 20 km w.e.. But a particular case of this uncertainty, namely difference between 
parametrizations for total cross section for absorbtion of a real photon by a nucleon at photonuclear interaction 
from Refs. J|J1, does not lead to a significant difference in resulting intensity. 



(c) Difference between muon spectra ( [2.1[ ) and (2.2) leads to uncertainty from -4% (D = 1 km) to 16% (D = 20 
km). 

(d) Error which appears due to simplified, entirely "continuous" ionization lies, commonly, at the level of 2-3 %. 
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(e) Dependence of simulated muon flux intensity upon v cut is the most weak one comparing with other studied error 
sources. Function I(v c ut) is almost a constant if v C ut < 0.02-0.05 and changes in a range ±1-2% which is very 
close to statistical error. 
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FIG. 5. Intensity of vertical atmospheric muon flux I at different dept hs D of pure water vs. v cu t as obtained by simulation 
with muons sampled according to sea-level spectrum from Ref. |B| (Eq. (2.1)). Closed circles: ionization is included in SEL; 
open circles: ionization is completely "continuous". Two horizontal solid lines on each plot show value for survival probability 
simulated with all muon cross sections multiplied by a factor k a = 1.01 (lower line) and k a = 0.99 (upper line) for v cut = 10~ 4 . 
Dashed lines on plots for D < 5 km corespond to intensity which was calculated for all energy loss treated as "continuous". 
Dash-dotted lines show intensity of vertical muon flux simulated with ionization included in SEL, v cut = 10 -4 and muons 
sampled according to Gaisser sea level spectrum (Ref. Eq. (2.2)). Horizontal dotted lines correspond to v cut — 10 -4 and 
cross section for absorbtion of a real photon at photonuclear interaction parametrized according to Ref. |}| (Eq. (2.4)) instead 
of parametrization proposed in Ref. fef (Eq. (f2.3|)) which is basic in the MUM code. 

We also tried to reveal how value of v cu t and model for ionization influence upon differential muon depth spectra. 
No differences were detected which would exceed statistical error. It is illustrated by two figures. Fig. [| displays 
simulated mean energies for vertical muon flux at different depths D of pure water vs. v cu t as simulated with muons 
sampled according to sea-level spectrum (f2~T). Two dependencies are presented at each depth for two model of 
ionization loss. No appropriatenesses are visible on the plot. In Fig. ^ we present simulated differential muon spectra 
at four depths D = 1 km, 5 km, 10 km, 20 km. Two spectra are displayed on each plot obtained i) for the most 
accurate value of v cu t = 10~ 4 and ii) for the most rough case v cut = 0.2. Spectra are normalized to 10 4 muons and 
at each depth are divided into three parts along muon energy to keep linear scale at Y-axis. Again, no statistically 
significant differences are marked. 
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FIG. 6. Mean energy for vertical muon flux at different depths D o f pure water vs. v cu t as obtained by simulation with 
muons sampled according to sea-level spectrum from Ref. || (Eq. (2.1)). Closed circles: ionization is included in SEL; open 
circles: ionization is completely "continuous". 
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FIG. 7. Simulated differential muon spectra at four depths D = 1 km, 5 km, 10 km, 20 km of pure water. Two spectra are 
displayed on each plot, namely simulated for v cu t = 10~ 4 (histogram) and v cu t = 0.2 (closed circles). Spectra are normalized 
for 10 4 muons with energies E > 0.16 GeV and at each depth are divided into three parts along muon energy to keep linear 
scale at Y-axis. 



9 




D (km w.e.) 

FIG. 8. Dependencies for relation I2/I1 vs. water equivalent depth in standard rock as computed in this work (closed 
squares), in Ref. || (open squares) and in Ref. [^] (closed circles). I\ is depth intensity for vertical muon flux simulated with 
ionization included in SEL, I2 is one sim ulated with entirely "continuous" ionization. Data for this work are obtained for 
sea-level spectrum from Ref. t| (Eq. (2.1)) and v C ut = 10_ -3 ; data from Ref. |^| represent result of simulation for spectrum 
from Ref. § (Eq. (|T§)) and v cut = 10" 
grouping" technique. 
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FIG. 9. Dependencies for relation I2/I1 vs. water equivalent depth in standard rock as computed in this work (closed 
squares) and in Ref. Q (open squares). 7i is depth intensity for vertical muon flux simulated with entirely "continuous" 
ionization and v C ut = 10 -3 , I2 is one simulated with the same treatment of ionization and v cu t = 10 -2 . Data for this work are 
obtained for sea-level spectrum from Ref. || (Eq. (2.1)); data from Ref. Ji) represent result of simulation for spectrum from 
Ref. § (Eq. 0). 



IV. DISCUSSION 



Results obtained in this work are evidence of accuracy in parametrizations for muon cross sections and spectra to 
be the principal source of uncertainties when simulating muon flux at depths where neutrino telescopes are located. 
Both for monoenergetic muons and muons sampled according to a power spectrum it gives uncertainty at least of 
2-4% to 10-15% in resulting intensity of muon flux, depending upon depth. Unfortunately, this level has at present 
to be considered as a limit for accuracy of muon propagation codes. Influence of model for ionization exceeds this 
limit only for monoenergetic muons with sea-level muon energies E < 10 TeV and only if level of observation is at 
very last stage of muon path where major fraction of initial muon energy has been lost. Actually, due to steep shape 
of atmospheric muon power spectrum, an essential part of muons reaches detector location being just on the last 
part of path. Therefore effect remains noticeable also for real atmospheric muons but in this case uncertainty was 
found in this work to be much less: 2-3%, which is in an excellent agreement with Refs. while Ref. [Ij predicts 

much more significant difference (up to 17%). We suppose this disagreement may result from the fact that "small 
transfer grouping" technique used for simulation in Ref. || treats muon cross sections to be constant between two 
interactions in contrast with algorithm used in the MUM code. As was shown in Ref. ||, in this case switching off 
the fluctuations in ionization leads to an additinal amplification in effective energy loss for muon energies E < 1 TeV 
and, consequently, simulated muon flux intensity must decrease at relatively small depths where muon spectrum is 
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formed by muons with sea-level energies just in a range E s ~ 10 2 -10 3 GeV. In Ref. || the same simplification was 
used but reported reults were obtained by simulation with v cu t = 10~ 3 . At this value of v cu t role of correct treatment 
for free path is not significant. Choice of value for v cu t is of even less importance and again, it is more critical if 
one investigates monoenergetic muon beam but for power spectrum alteration in v cu t within v cu t < 0.02-0.05 leads 
only to 1-2% differences in simulated muon flux intensities which, again, is in a good agreement with level of errors 
reported in Ref. || . Differences between muon flux intensities simulated for different models of ionization and values 
of Vcut, as obtained in given work and in Refs. are presented in Fig. || and Fig. |[ 

Since computation time which is neccessary for simulation depends strongly upon used model of ionization loss and 
upon v cut , it seems to be reasonable for most purposes to set v cu t to v cu t ~ 0.01-0.05. Anyway, in this case resulting 
error will be less comparing with error caused by accuracy of parametrization for muon cross sections and sea-level 
muon spectrum, if one treats free path L by an accurate method (Ref. || ) . Note that "inner" accuracy of the MUM 
code was found to be worse than "outer" one (which results from uncertainties of muon cross sections) if ionization 
is set to "continuous" model and v cu t > 10 -2 . For fluctuated ionization inner accuracy remains high enough at least 
till v cut = 0.05 (Ref. So, using MUM one should set v cu t to v cu t < 0.01 if ionization is entirely continuous and 
to Vcut < 0.05 for ionization included in SEL. For test runs it is possible to set v cu t even to larger values. On the 
other hand for some methodical purposes it may be useful to simulate fluctuations in knock-on electron production 
and choose more fine v cu t, e.g., if one wants to exclude an additional error when comparing results of simulations for 
different models of muon sea-level spectrum with each other. 

It is impossible to consider all particular cases and give some strict conformity between setting of parameters at 
muon MC propagation code and problem to be solved. We had for an object to create a code which would allow to 
change easy the most important parameters and to perform an analysis which would allow anyone to choose these 
parameters according to one's purposes and taste. We hope this object has been achieved with this article. 

V. CONCLUSIONS 

We have presented the detailed investigation for dependence of computed depth muon flux upon the most important 
parameters of muon propagation MC algorithm, which was done with the MUM (MUons+Medium) code (Ref. 0). 
Contributions of different simplifications and uncertainties to the resulting error have been analysed and ranked. We 
hope that presented work can be useful for further development of muon transport algorithms which are neccessary 
for adequate analysis of muon data obtained at underground/water/ice neutrino detectors. 
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